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INTRODUCTION 

The following sections of this final report on high 
speed hydrogen/graphite interaction deal with the various 
aspects of the research program. 

Section 1 concerns the theoretical predictions of 

hydrogen/graphite erosion rates. Our previous report* 1 ' 1 ) 

covered this topic area in some detail, however, this phase 

of the research is continuing under Dr. Sharma's direction 

and a brief summary of the activities are included in this 
section. 

Section 2 describes high temperature, nonequilibrium 
hydrogen flow in a nozzle, in this section some detail as 
to computer programs and subroutines is supplied in order to 
calculate flow characteristics of nuclear rocket applications. 
The analysis includes subsonic, transonic and supersonic 

regimes together with equilibrium and non- equilibrium con- 
siderations. 

Section 3 describes the molecular beam studies of 
hydrogen/graphite erosion. This phase of the research has 
concluded. The findings, together with a description of the 
approach and experimental apparatus. 


are summarized here. 


liijiHiijjiJijj'aiaiil 
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SECTION 1 

— m p B .nTPTTQ HS OF HyflP OPEN /GRAPHITE EROSION RAT ES 

A theoretical model ie developed to predict the erosion 
rate of graphite coating by a stream of hydrogen gas. It 
based on a phenomenological description of turbulent reacting 
boundary layer flow. It also takes into account the effect of 
lateral curvature in the nozzle as well as the dissociation 

of hydrogen gas. 

The basic conservation equations with the relevant bound- 
ary conditions for this model were summarized in the last re- 
port' < 1 - 1) . The formulation of the problem was done by consider- 
ing the multicomponent diffusion of species. The detailed de- 
scription of the transport properties needed for the solution 
of the equations as well as the transformed forms of^the gov- 
erning equations were also given in the last report 

A numerical solution is being sought for the set of 
coupled nonlinear partial differential equations. For this 
purpose, we are following the general procedure adopted by 
Herring and Mellor for calculating the boundary layer flow 
of a single component (non- reacting) mixture. Using th 
Michelson scheme, the partial differential equations have been 
transformed into finite difference equations which are then 
being solved by Gaussian elimination procedure. The expression. 

are too lengthy for presentation in this report. 

Thete ntially three successive Phases invol 

solution of this problem: 


the complete 
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(a) The flow external to the boundary layer must be 
known. This task was undertaken by Mr. R. Hamman who provided 
the numerical data for the non-equilibrium flow field with a 
given nozzle geometry by neglecting the boundary layer effects 
(see Section 2) . The choice of step- size along the axial direc- 
tion is consequently restricted by this input data. 

(b) In order to solve the governing equations step by step 
along the axial direction, we use approximate turbulent similarity 
flow equations at the initial station. Herring and Mellor gener- 
alized the results obtained from incompressible equilibrium tur- 
bulent boundary layer flow studies to get approximate compres- 
sible turbulent similarity flow equations. In the present anal- 
ysis, we make use of the same prescriptions for the momentum 
conservation equation expressed in terms of temperature in con- 
trast to enthalpy. Moreover, the same concept is extended to 

the species conservation equation. 

It must be pointed out that since we have included the 
dependences on temperature and composition for the various 
physical properties of the mixture, the final expressions are 
far more complex, and the logic of the computer program also 
becomes very involved. 

(c) Having obtained the results for the similarity flow, 
the calculations can be easily completed for other stations 
by considering the general form of the governing equations. 

No essential difficulty can be expected at this stage. All 
the relevant details have been worked out. 

In view of the unexpected change of location of the in- 


\ 
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vestigator for this phase of tne program. Dr. 0. P. Sharma, 
and the necessity of using a smaller computer facility 
(IBM 7044) , the debugging of the Fortran has been rather 
slow. It is anticipated that the computation will be com- 
pleted by summer 1974 at no additional expense to the contract. 
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SECTION 2 

HTftH TEMPERATURE. NONEOUILIBRIUM HYDROGEN FLOW IN A NOZZLE 
2.1 Analysts 

The analysis is composed of three rather distinctive flow cases: subsonic 

one dimensional flow, transonic flow and supersonic axisymmetric flow. 

Only a rough outline of the respective theories is given. For more de- 
tailed information one is referred to the references for this section. 

The subsonic flow analysis consists of a pressure defined one-dimensional 
non-equilibrium analysis as given in Refs. 2-1 S 2-2. Insufficient reliable 
information was known to warrant a computer program on a two-dimensional 
axisymmetric treatment of the subsonic flow. 

In the case of transonic flow a model is chosen for a simple, more geo- 
metrical treatment analogue to that discussed in Reference 2-2 which gives 
accurate results compared to experiments. An extensive treatment of the 
transonic axisymmetric case is given in Reference 2-2, but was rejected 
for this program as being too complicated. 

The analysis of the axisymmetric non-equilibrium supersonic flow employs 
a method of characteristics solution for the partial differential equations. 
The finite difference scheme chosen stems from Reference 2-4, being a 
more simple method than that used in Reference 2-5. 

2.1.1 Subsonic one -dimensional flow 


The governing differential equations for one-dimensional flow with chemical 
non-equilibrium are 




sw 


(tewtiwvittj ^ 
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In the disoussion we will not consider the frozen case, as its method 
of solution is almost identical to the one for equilibrium flow. 

As the set of equations both in the equilibrium and in the non-equili- 
brium oase are to be integrated through the transonic region, the area 
ratio oannot be used as the independant variable because the area ratio 
defined set of equations is singular at the sonic point ' 

Using the velooity or pressure as the independant variable has the extra 
advantage that the set of four differential equations, which originally 
required simultaneous integration, is reduced to a set of three equations. 
_With the flow field the corresponding area ratio can be computed afterwards. 


In the non-equilibrium case there exists a coupling of the equations with the 
area ratio by means of the chemical relaxation equation. This case will 
hare to be solved by iteration of the complete subsonic flowfield until 

the computed area ratio converges to. the area ratio defined by the nozzle 

shape. 
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2. 1.1.1. Equilibrium flov 


Using velocity as the Independant variable equations (2.2) and (2.3) 

4 * 

can tie written 


5 W dv el 


p * © 





4a.6) 
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Equation (2.4) is replaced by the condition for equilibrium or- w^(p,T). 
As the enthalpy h is a function of o£ and T only 


where 


c4U« U^.dT 




ORIGINAL PAGE is> 
OF POOR QUALITY' 


Combining this result with (2.5) 



4km k T olT It * £ clp -v *‘^2±a dT J 

1 •*T’ 


which reduces the set of equations to 

« , 

£W«Ut -e- alp m o 


(2»ft) 


Writing this in finite difference fora with the Telocity as the independant 

variable 


/ 


Ap ■ £ _jw ) Aw 






which combined with the condition for chemical equilibrium 


* r *cCp/0 


gives a solution for the equilibrium flow properties. 


The coefficients in (2*9) end (2*10) are taken to be the average values 
between two consecutive points, while the density is computed using 
the equation of state 


3 2» h 




The local frozen speed of sound can be found using the equation given in 
appendix A. 

As the mass flow m is constant, it follows from 

f Aw* N 

that at the throat, where A - A*- A nin , jw should take a maximum value. 
After location of this maximum, the throat properties and the equilibrium 
aasa flow can be determined, and the area corresponding with every set of 
properties p, T, and w is known: 

*h* a% u 



An sxpression for the nozzle shape links the area with the 
which permits us to use the functions 


U-w.) 

axial coordinate, 
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w» (,-«-) 

T, T h i-K) 

• fit«^ f-) 

and all derived equilibrium properties, whenever they are needed. 

2.1. 1*2. yon-equilibrium flow 

As the -pressure is relatively insensitive for non-equilibrium effects, it 
is convenient to use it as the independant variable in the kinetic analysis. 

Equations (2.2), (2.4) and (2.3) can be written 


l^wjolw + [1? j dx. =■ o 
dot s o 

or 

* ,w -' C 

dec w dhc 

dT. f- b « V _L 

The pressure and its derivatives are assumed to be known; as a first ap- 
proximation is taken to the equilibrium pressure distribution. 

A fourth equation can be derived by differentiating the equation of state 

(2.11) 


«i*. [iirjrj. 14T . _l 

*• PValX.^ T* doc »-r«c J 


These four equations can be rewritten, an shown in appendix A, in the form 
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where 


, f - =i l4»\ 

dot ** y-ci'H.J 

l2-»0 

£*S » K r ii 

CiX. * Vi/ 

(?•*) 

. r iii (.£?). b i t 

oi>w L yp v cix./ J 

l*-*0 


C 5 '7) 

• 

b« ia t. 

-n. T 11 

• 

A, it _b 

ft 


l+»et 


The integration method ia taken from Reference 2-1 with a few minor 


alterations. 


Setting 

W 

Sk'* 

V* 

u, » 

Act 

•s* 4? 


and 







four second order integration formulae in the increments of Y^ can he de- 


i 


rived of the form 
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where 


M , 

4^1 ^V»js* ^\».*** , *’« ) ~- H - ( 4-W h ) 1 

U *»r« J 


ll l m etepsize for the 4 th interval. 

The function value, f . and the derivatives * and £. . are taken 
_ th 1 * J 


th . A 
n point. 


at the 


4 * t ’** , *“ inc ™'"‘ *•»•« th. computed increment h.e to no knevn te 
»pply equation (2.18), it i. nec.re.ry to nee a different act of formula, 
for the starting step of the integration. 

The first order variant of equation (2.18) is 

Ki % *>. k £ A . . L. \ U /., 

'* +»• ix, |0 vi*o k j,i ; K 

when applied on the first interval. 

The relative error in th. increment , or ti„ e.cond order method i. 

A Ki» I -I,..., I 

1 ale. -lc. I 

Ihi. fonmil. i. used to determine the allowable atepeise. 

A. error v.ri.a a. the .fp.it. aquared. the criterion for donhlin, 

•nd halving the steps ize is given as 


L* 

»o 

< < s 

(,2.20 a*) 

0 l 


At H. >S 

U-2«kJ 

. o L 
* %*%«•! 


A W • ✓ s* 

to 

[z.ioc ) 


where & is the maximum 


error allowed. 


Equations (2.18) and (2.19) can he written as four linear algebraic equa- 
tions 

M t v -, *Vt ♦•it *«-i 

e # 

which are solved for the increments k. , , i»l,2,3,'». 

i,n+l* 1 11 

The integration method is only weakly dependant on the engine sise and 
allows an integration step several orders of magnitude larger than the 
chemical relaxation length. 

The gas properties at the throat and the mass flow are determined the same 
way as treated in the paragraph on equilibrium flow. 

♦ 

** * A * 

and 

A u. 

A convergence criterion is applied on the mass flow and the one— dimensional 
computation is ended, when sufficient convergence is achieved. If no suf- 
ficient convergence is found, the complete one-d imonsional non-equilibrium 
calculation is restarted, using the newly computed pressure profile. 

The form of the functions and derivatives in equations (2.18) and (2.19) 
is given in Appendix A. 

Transonic flow 


Basis for the analysis is a modified Sauer analysis as given in Reference 22 
by Nickerson. The method employs first order perturbation theory to find 
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* solution for a constant property line at the throat minimum point. 

® more extensive discussion of the theory us€ Re ference 2—2. 
Only the results of and changes in that theory will he given here. 


The constant property line selected is the isohar at the throat minimum 
point. The coordinates of the isobar are given bv 


where 



c 





a* ) Vi 

r+i 


(?•*') 


Uj - downstream normalized throat radius of curvature 
Y - average specific heat ratio in the throat region 
r ■ radial normalized coordinate 
* • axial normalized coordinate. 


The average specific heat ratio is determined as 



ii.Zi.) 


where the asterisk denotes the throat property following from the one-dimen- 
sional analysis and is the property at point i on the isobar, n is 
the total number of points on the isobar. 


The spacing of the points on the isobar is determined such that the mass 
flow through each annular element is equal. 

Hence 


and 



(a.«; 

(/•*».) 



-**•*»• — (u.l) (dt.zrj 
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Denoting all one-dimens innal throat properties with an asterisk, it follows 
for the pressure on the isobar 


I L. ) 


The axial and radial components of the velocity are 




^ ** V-t-i -J ' 


v * ^ l J l ^ i 


Using equation (2.23) these expressions can he simplified into 


Uta* (• *+ JL \ 

1 liC. / 


>->■>[ 


which shows that the axial velocity is constant on the isobar. 


The resultant velocity is 


UWJ (2.11 

The streamline angle is found easily by 

e.atfcm (.■£) (2.3C 

The two other flow properties to be determined are the mass fraction atomic 
hydrogen and the temperature. 

As solving the chemical relaxation equation 


4e( u» 
Sx " W 


over the relatively large distance from throat to isobar leads to 


severe 
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numerical instability for interpolation from the property table provided 
by the one-dimensional non-equilibrium analysis. The temperature follows 


from the energy equation 


u. -u l 


Description of program subroutines 


(*•*0 


‘ xhe program subroutines are divided in four groups 

- main program and general subroutines used throughout 
the computation 

- subroutines used for the computation of the one-dimensional 
subsonic flow 

- subroutines for the determination of the transonic line 

- subroutines that perform the two-dimensional supersonic 
flow calculation. 



| 

\ 


The subroutines in this group are 

i 

i 


MAIN PROGRAM 
SUBROUTINE C0NSTA 
CHEMEQ 
EQUIL 
GAMMA 
TEMP 



MAIN PROGRAM 


<■ o - 


4.2 Input and physical constants 

SUBROUTINE C0NSTA 

In this subroutine the physical constants used in the program are 
initialized and the nozzle shape and case parameters are read in. 
The variables, used in the internal computation, are derived 
from the data given in the nozzle shape parameters. 

Finally the values of the input parameters are printed. 

Physical constants 


The physical constants, as they are given in the subroutine, are 

- AA constant in the exprsaion for the reaction rate PKF- 

, where T is the temperature in degrees K. 

- fMK mass of a hydrogen atom 

- FK Boltzmann constant 

- TOED characteristic temperature for dissociation of 
hydrogen 

- TOEV characteristic temperature for vibration of hy- 
drogen 

- PI geometrical constant 

- WHY atomic weight of hydrogen 


The constants are set at the values 

AA ■ 0.7x10*® cm^ °K^mo lessee"* 
• fMH - 0.l674xl0“ 25 g 

FK - 0.138xl0" 17 erg °K~ 1 aole“ 1 

TOED • 54733 °K 

TOEV - 6140.0 °K 

PI . 3.14139 

WHY • 1.008 g/raole 


A 16 - 


— IF ( LENGTH. GT.O) it follows through 

RCHIT- • REE -TMEX.ton^THJUp) 
A1 x — 

TMEX* 

* *t«wk (TMSOF*> J 2 . A*1 . XCC 


- XEP » XC.fi 4TH6X 

• TM6X * otOM (l.AH.xEP) -*-B1 (rod) 


- IF ( LENGTH. l^.Q) 


2 ( fifix » r. i _ cTHfifi. (-1 .cos *})) 

• XEFt L rXEE 

ton (Thex) 


. COL » 


(xer -xfifi)towtTHtuF ) 

(THSUP*b .to** irvtEx) 


The parameters I WALL, IPEKF, IPROP, INTER, NPNTS, IOUT and LENGTH deter- 


mine the case that has to be computed and the shape and quantity of the output. 
PROGltAM INPUT Pfi stagnation pressure (dyne/cm^) 

T fi stagnation temperature (°K) 


PROGRAM OUTPUT 


internal constants 
PRINT 


CALLING 


DATA 


SUBROUTINE EyilIL 


CALLED BY 


General subroutines 


MAIN PROGRAM 

SUBROUTINE CHEMEQ 
EQUIL 

GAMMA 


TEMP 
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Nozzle shape parameter a 

The nozzle shape ie assumed to be composed out of circular arcs, straight 
lines and a parabolic supersonic expans ipn section, 

I 

Head in are 

• KC1IAM combustion chamber radius 

- K0NE nozzle inlet radius of curvature 

— KTWjf upstream throat radius of curvature 

- KT1DIE downstream throat radius of curvature 

- KSTAR throat radius (cm) 

- 1USXIT nozzle exit radius 

— TIISUB maximum nozzle inlet angle (rad) 

- THSUF maximum nozzle exit angle (rad) 

- THEX IP (l£NG'Ai.Ey.O) nozzle exit angle (rad) 

IF (LENGTH. GT.O) length of parabolic nozzle section 

All values, with the exception of KSTAlt and the angles, are dimensionless 
with respect to the throat radius. 

The internal nozzle shape parameters, which can he derived 


X&E m NC . A»V* 

- KCS. r RGM*M --1 - ATwtf (l - C.O V HSuO,')) 4- 

1 

+ - 1 > 

t«v%(TMSUO) 

• KDE. xce -t- «cr wjff . aiM (TH^va) 

- X EE» XOE ♦> RTHftE . AIM 

• GEE . RGMAM * (tot -"I ) 

.Rt£ . -1-*- RTW0 (1-COS (THSUB) ) 

- REE. * -d * ttT HRE (i- 
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SUBROUTINE Cll»ir:u (T.FKP.DFKP) 

This subroutine is used during the complete computation to provide the 
other subroutines with the value of the equilibrium constant FKP ((dyne/ 
cm 2 )^) and the derivative 1)FKP (-—(ln(FKP))) at a given temperature T 
for the reaction Ilg!* 131. 

The first time the subroutine is called two arrays, S and TAB, are read 
in. Immediately afterwards a flag is set 
NFLAO-1 

which prevents that the next time, when the subroutine is called, it will 
try to read the values for S and TAB again. TAB and S are defined 

- TAB ■ ***log Kp , Kp in (atmospheres)^ 

- S corresponding temperature in degrees K. 

and are taken from the JANNAF tables, Reference 2-6. The range of tempe- 
ratures covered runs from 0°K until 6000°K. 

At the required value of the temperature T the corresponding value of FKP 
is found by logarithmic interpolation, and the result is converted to the 


t 

output units. In 
determined. 

the same 

loop the value of the derivative — =r(ln K ) is 

via r 


PROGRAM INPUT 

T 

temperature (°K) 

1. 

PROGRAM OUTPUT 

FKP 

equilibrium constant ((dyne/cra 2 )^) 

1 

i 

t 

t 

CALLING 

DHCP 

logarithmic derivative of the equilibrium 
constant ((°K)”*) 

DATA 

* 

CALLED BY 

SUBROUTINE 


i 

j 



< 
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SUMtOUTINK KUUIL (P.T.ALP) 

Calculates at given pressure P and temperature T the corresponding value 
of the mass fraction atomic hydrogen ALP at chemical equilibrium con- 

ft 

ditions according to the relation 

ft g , 

ALP-DSqUT( 1 ./( 1 .+4 . *P/ ( PKPnFKP) ) ) 
where FKP is the value of the equilibrium constant at T °K. 

PKOOHAM INPUT P pressure (dyne/cm 2 ) 

I temperature ( °K) 

PHOGltAM OUTPUT ALP mass fraction atomic hydrogen at equili- 

brium conditions* 

CALLIN G 
CALLUP HY 


sumjAJTiNK cinMtq 
SUIR^UTINE 


SUatdUTINE GAMMA (P.T.ALP.CAM.AF2.KHiO 


At given 

- pressure P 

. - temperature T 

- mass fraction atomic hydrogen ALP 

are computed 

- specific heat ratio GAM of the mixture 

- squared frozen speed of sound AK2 


- density Wtyf 


according to the 

expressions 

(Al) , (A7) and (AH) of appendix A* 

PKOGHAM INPUT 

P 

pressure (dyne/ cm**) 


T 

temperature (°K) 


ALP 

mass fraction atomic hydrogen 

PK OGHAM OUTPUT 

GAM 

specific heat ratio 


AF2 

squared frozen speed of sound (cm^/sec^) 


Wtyf 

*1 

density (g/cnr) 

CALLK1) BY 

SUHK0UTINE 
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2.5 Subsonic flov 

The subroutines in this group are split in two subsections, the first 

» 

one dealing with one-dimensional equilibrium flow, the second one with 
one-dimensional non-equilibrium flow. It contains 

| SUHH/fUTINK C/to} 

CffaTAJ* 

ST0UT1 

1 TABLE 

CHKCK 

PERISH 

PUN/fUT 

SffLVE 

SVliCffl 

SUBC/^2 

SUBKX 

Equilibrium flow , 

SUBROUTINE C0EQ 
C0RTAB 
EQUIN 
ST0UT1 
TABLE 


t 

~/A » w w r- 
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SUHRglJTINE CtfeO (P.T.W.ALP.C01.C02) 

The coefficients of the two finite difference equations for integration 
of the one-dimensional equilibrium flow field (c.f. section 2. 1.1.1.; 
equations (2.9) and (2.10)) 


CO?4 > 




I 

V -gg .? - \ 

L ^ , 


v 


are computed here. The expressions for the partial derivatives 



*vr 


can be found in appendix A (equations (A12a) and (A12b)). 


PHOOIUM INPUT 


FHOGKAM OUTPUT 


CALLING 


P pressure (dyne/cm 2 ) 

T temperature ( °K) 

W velocity (cm/sec) 

ALP mass fraction atomic hydrogen 

Cql coefficient in equation (2.9) 

coefficient in equation (2.10) 

SUHtjfclTINE ClQMEq 


CALLUP BY 


SUHH0UT1NL E(JUIN 
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SUHHrfUTLNE C^KTAB (W.X.M.WK.KX.J) 


The eubroutine modifies a table (W,X) with M entries into a table (WE, 
EX) with J entries with the result that hetween each, two values of the 
array EX in the second table a constant interval is constructed. 


The interpolation between two values of X is achieved with the formula 


weta&r W(X) CL 


t ^ w \ 

K H t 


X (t) <EX < X^X 


where EX(j)«EX(j-l)+UELL 

RELU.03*KSTAlt 
EX(1)-X(1) 


A logarithmic interpolation method is chosen as it represents most ac- 
curately the shape of the function W«W(X) between two subsequent points. 


PROGRAM INPUT 


PROGRAM OUTPUT 


CALLED BY 


W velocity (cm/sec) 

X axial coordinate (cm) 

M number of entries in (W,X) table 

WE velocity (cm/sec) 

EX axial coordinate (cm) 

J number of entries in (KE,EX) table 

SUBROUTINE E^UIN 
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SUBROUTINE KljUIN (Pp.TlUMASS.M.X.P.re.WS.AB.PlEQl 

The subroutine computes the equilibriua pressure distribution p as a func- 

* 

tion of the axial coordinate X and the equilibriua mass flow MASS. 

Two nain calculations can be distinguished* After the computation of the 
initial flow conditions from the stagnation pressure P p and stagnation 
temperature T 0, which uses the comnon expression for the temperature ratio 
in compressible flow as a function of the Mach number and the specific 
beat ratio, the energy equation and the condition for chemical equilibrium, 
a table of the equilibrium gas properties 

- velocity V 

- pressure P 

- temperature T 

- mass fraction atomic hydrogen ALP 

- density WI0 

- frozen Mach number HUGH 

is computed, using the relations described in section 2.1. 1.1. 

During thit calculation the maximum value of the product (RHflmV) is de- 
termined, from which the equilibrium mass flow MASS follows. 

After computation of the sets of gas properties corresponding with the 

area by applying continuity, follows the construction of a table of the 
velocity WE as a function of the axial coordinate EX by means of SUBROU- 
TINE TABLE, which employs an analytic expression for the nozzle geometry. 

At the same time the nozzle inlet conditions 

- pressure PBKG-PS 

- temperature TS 

- mass fraction atomic hydrogen AS 

- velocity WS 
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SUBRrtJTINK SlgUTl (P.T.ALP.iaK>.W.A.X.^Ant) 

Thia subroutine provide# the regular output of the result® of the one- 
dimenaional performance calculation. 

The nozzle area is converted in the area ratio by dividing by (RSTAK* 
KSTAKaPl) . Printed are for all stations 

area ratio A 

■ axial coordinate X (cm) 

- pressure P (dyne/cm 2 ) 

- temperature T ( °K) 

- density Hltjrf (g/cnr 5 ) 

- mass fraction atomic hydrogen ALP 

- velocity W (cm/sec) 

~ frozen Mach number FMACH 

When the frozen Mach number becomes greater than unity, the vacuum specific 
iapulie it computed according to 




and printed. 



PHOGUAM INPUT 

P 

pressure (dyne/cra 2 ) 

• 

7 

temperature (°K) 


ALP 

■ass fraction atomic hydrogen 


HH0 

density (g/cm^) 


W 

velocity (cm/sec) 


A 

nozzle area (era 2 ) 


X 

axial coordinate (cm) 


FMACH 

frozen Mach number 

PROGRAM OUTPUT 

PRINT 


CALLED BY 

SUBROUTINE EQUIN 


SUBEX 
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i 

i 
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i 

i 


i 

i 

f 
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The table (WK,EX) is corrected with the ai.d of SUBROUTINE C^llTAB to. en- 
sure an equal stepsize DX in the resulting table (W,X), which furnishes 
the new input for the second run through the equilibrium flow integration. 
The final pressure table (l\X) is obtained again with the aid of SUHR0U- 

4 

TINK TABLE. 

SUBROUTINE ST^UTl takes care of the standard print output of the results. 


PROGRAM INPUT 


PROGRAM OUTPUT 


CALLING 


P 1 
T ft 

niAas 

rt 

X 

p 

TS,WS,AS 


stagnation pressure (dyne/cm 2 ) 
stagnation temperature (°K) 

equilibrium mass flow (g/sec) 
number of entries in the arrays P and X 
axial coordinate (cm) 
pressure (dyne/cm 2 ) 

temperature (°K), velocity (cm/sec) and mass 
fraction atomic hydrogen corresponding with 
the first entry in the (P,X) table 


PBKG nozzle inlet pressure (dyne/cm 2 ) 

SUBROUTINE CHTMEQ 

cOw 

/ cOltTAB 
B1UIL 
GAMMA 
STOUTl 
TABLE 


CALLED BY 


subroutine suhex 
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SUlttidUTlNK TABUS (P.A.N t AST A K.T.W.AlP t X t g t M t TS,V , S 1 AS,NSTAK j - 

The property/area table (P,A) with N entries ie converted in a property/ 
axial coordinate table (q,X) with M entries. 

The values of the other three one-dimensional flow properties IS, US 
and AS corresponding with the first entry in the (q,X) table are inter- 
polated from the input arrays T, V and ALP. NS TAR is the index of the 
entries at the nozzle throat. 

The relations that link the axial coordinate X with the nozzle area A, 
non-dimens i ona 1 i z ing the axial coordinate X and the radial coordinate 
by dividing through the throat radius USTAR, are: 

k(J)«VMJ)/astak 1 

- if (R(J).LT.RCHAM) 

X( J)« ^(RfJNlS 2 -(U( J) -RCHAM+H^E) 2 ) 

- if (R(j).LT.RBi:) 

X( J)-(ltCHAM-U( J)+K0NE( cqs (figuiiy tanCmsUBT 

- if (H(J).LT.RCE) 

X(j)«XDE -^(UTW0 2 -(UTW/+ 1-R( J ) ) 2 ) 

- if (J.OT.NSTAR) 

X(J)-XDE +^(R'niRE 2 -(R(J)-l-RTia« l 2) 2 ) 

- if (r(j).gt.rke) 

x( J)«XKB ♦ C0L( 1- \f( } 

Finally all computed X values are dimens ionalized 
X(j)«X( J) .USTAR 


fftOGltAM INPUT 

PA 

* gai property/nozzle area table to be con 

, a 1 ^ 

Ter ted; A in cm 


N 

number of entriee in (F f A) table 

o 
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N entries 


NSTAlt 

index of properties at. the nozzle throat 

HtOGUAM OUTPUT 

y.x 

gas property/axiftl coordinate table; X 
in cm 

* 

M 

number of entries in (q,X) table 


TS.WS.AS 

with (q(l),X(l)) corresponding values of 
the arrays T, W and Ali* 


CALLED HY SUJIifAJTlNK EqUIN 

SUHEX 


3*2*2« Non-equilibrium flow 

SUDUjfUTINK CHECK 
DBtlSH 
PUNjfUT 
S0LVF, 
SUHC/Jl 
SUBC02 


StJBEX 
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SUBROUTINE PERISH ( P ,T .RHO.ALP.W . DPX . DP2X . K1 . FI> . F3 1 F4 . FIX . FI V . FIR. FLV. 
Fat . P2T .F2A.F3X.FW. F 3 H . F 3 T . HA . MX A FW . F'iH . F4 T . Ft A) 


This subroutine computes at known 

- pressure P 

• - temperature T 

- density lQlfi 

- mass fraction atomic hydrogen ALP 

- velocity W 

- pressure first derivative DPX 

- pressure second derivative DP2X 

the functions FI, F2, P3, F4 and their derivatives with respect to x, 

, T and , according to the expressions given in appendix A (eq- 
uations (A13) through (Al6) and following), which will be used for the 
calculation of the coefficients of the first and second order integra- 
tion methods for one-dimensional non-equilibrium flow. 


PROGRAM INPUT 


program oumrr 

CALLING 
CALLED BY 


p 

pressure (dyne/cra^) 

T 

temperature (°K) 

uhO 

density (g/cur*) 

ALP 

mass fraction atomic hydrogen 

V 

velocity 

DPX 

pressure first derivative (dyne/cm^) 

DP2X 

pressure second derivative (dyne/cm*) 


F1,F2,F3,F4 and their derivatives (c.f, appendix A) 

SUBROUTINE CUttlLV 

SUBROUTINE SUHCOl 
SUBCff2 


1 


- 30 - 

SUBK^UTINli SdLVB ( P . f .Ittlrf . ALP . W . 1)X . DUX . DP2X . JftlN , PTN . 1>HN . PAN . PV .PT.DH, 
DA.NdNK.PXl) 

Four linear algebraic equations in. four unknowns are solved using Kramer's 
itule. The equations are of the form 

a^W ♦ bjI>A ♦ ^PT ♦ djPH - e t i - 1,2, 3, * 

and are the rearranged integration formulae from section 2. 1.1. 2, equa- 
1 tiona (2.18) and (2,19). 

The coefficients a^, b^, c^, d^ and e^ are computed in SU19tyft!TINK SUBC01 
(first order method, applied only for the initial integration interval, 
if (N0NE.EQ.1)) and SUBUj^UTINE SUHCf$2 (second order method, used for 


all subsequent steps). 

• 


PKOG'.iAM INPUT P 

pressure at last calculated (n^* 1 ) point 
(dyne/cm“) 

T 

temperature at n** 1 point (°k) 

Hlfjtf 

density at n** 1 point (g/cm^) 

ALP 

mass fraction atomic hydrogen at n** 1 point 

V 

velocity at n** 1 point (cm/sec) 

nx 

axial distance between n and (n+l) 


' point (cm) 

DPX 

value of ^jjj* at the n** 1 point (dyne/cm^) 
value of at n**’ point (dyne/cra* 1 ) 

DP2X 

DWN 

dx“ th 

increment in velocity between (n-l) and 


n** 1 point (cm/sec) 

s 

DTN 

increment in temperature between (n-l)*’*’ 
. and n** 1 point (°K) 

DHN 

increment in density between (n-l)** 1 and 
n** 1 point (g/cm^) 


PROGRAM OUTPUT 


CALLING 


- 31 - 

RAN increment in mass fraction atomic hydro- 

gen between (n-l)**’ and n* h point 
N^NE n 

J1X1 axial distance between (n-l) th and n th 

point (cm) 

IIW increment in velocity between n** 1 and (n+l)** 1 

. point (cm/sec) 

DT increment in temperature between n** 1 and 

* 

(n+l)** 1 point (°K) 

DR increment in density between n and (n+l) 

point (g/cm^) 

PA increment in mass fraction atomic hydrogen 

between n** 1 and (n+l)** 1 point 

SUBROUTINE SU1JC01 
81JHC02 

SUBROUTINE SUBKX 


CALLED BY 
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SUHKdUTINK SUKCd 1 ( P ,T .KHtf .ALP . W . DPX . D?:ix . A l . A2 . AT . A4 . »;* . m . H’» . C2 . C3 . 
C*.D1.D2.D3.1>4.E1.E2.E3.E4.1«) 

Computation of the coefficients of the four linenr Algebraic equations, 
obtained from the first order integration method, c.f. section 2, 1.1. 2, 
equation (2.19)* 

The values and the derivatives of the functions FI, F2, F3 and F4 at 
given P t T, HH j6, ALP and V (equations ( A13) through (A18), appendix a) 
are computed in SUHH0UTINE PEltlSB. 

The coefficients Al, A2, etc. Are defined 

A1-F1V-1/DX : ' 

A2-FJV 

A3-F3W 

A4-F4V 

B2-F2A-1/DX 

BJ-F3A 

B4-F4A 

C2.F2T 

C>F3T-1/I)X 

C4-F4T 

D1-F1K 

D2-FJSI 

D3-F3K 

D4-F4H-1/DX 

EN-F1-F1X.DX 

EShi»F2 

K3.-F3-F3X.DX 

E4.-F4-F4X.DX 


< 



PROGRAM INPUT 

P 
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pressure (dyne/ cm") 


T 

temperature (°K) 

• 

RUf* 

density (g/cnr*) 


ALP 

mass fraction atomic hydrogen 


V 

velocity (cm/sec) 

4 • 

DPX 

derivative dp/dx (dyne/cm 5 ) 


DP2X 

derivative d 2 p/dx 2 (dyne/cm 4 ) 


DX 

axial interval (cm) 

PROGRAM OUTPUT 

Al ,A2, « 

••,K3»E4 coefficients of linear algebra 



equations 

Galling 

suretjfuTiNK deuish 

CALLED HY 

SURtt)*UTlNK SjtfLVE 


C» .1)1 . D2 . D3 . D't . E 1 . K2 . K3 . Kk . 1)X . 1>WN . 1)TN . PAX . D1LV . PX 1 ] 


Computation of the four linear algebraic equations, obtained from the 
second order integration method, c.f. section 2. 1.1.2, equation (2*18). 

jtht values and' the derivatives of the functions FI, F2, F3 and F4 (e- 
, quations (A13) through (A18), appendix A) are calculated in SUHU^UTINE 
^ DEKISB. 

Setting 

p«dxi(px+pxi)/dx 

the coefficients Al, A2, etc. are defined 

Al-FlW-l/Fq 
B2-F2A-1/PV 
C3-F3T-1/P<i 
D4-F4K-1/PV 
El— UWN/q-FlX .PX-F1 
Ea.-DAN/q-F2 
E3— DTN/«i-F3X . PX-F3 
Et— DUN/q-F4X. DX-F4 

| A2, A3» A4, H3, B*t,‘C2, Ck, Dl, D2 and D3 have the same definition as 

j In SUHiffUTINE SUHCjfl. 

| PKOGHAM 1KPUT P pressure (dyne/cm^) 

| , T temperature (°K) 

1 Itftyf density (g/cm^) 

- ALP mass fraction atomic hydrogen 

W velocity (cm/sec) 

• PPX derivative dp/dx (dyne/cm^) 



nx+nxi] 







PROGRAM OUTPU T 


CALLING 
CALLED HY 
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v 


DX 

DWN 

DTN 

DUN 

DAN 

DX1 

A1,A2 


Axial interval between last computed point 
and next following point (cm) 
increment in velocity between last computed 
point and the point before (cm/sec) 
increment in temperature between last com- 
puted point and the point before (°K) 
increment in density between last computed 
point and the point before (g/cnr) 
increment in mass fraction Atomic hydro- 
gen between last computed point and the 
point before 

axial interval before DX (cm) 

coefficients of linear algebraic 
equations 


SUBROUTINE PERISH 
SUBROUTINE S/JLVK 


Transonic flow 

This group is composed of only one subroutine 


SUBR0UTINE INLINE 


I 
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SECTION 3 

MOLECULAR BEAM STUDIES OF HYDROGEN/GRAPHITE EROSION 
Introduction 

As part of its contoured exhaust nozzle, the NERVA solid 
core nuclear rocket flight engine (Ref.3-l)was designed tc incorporat 
a radiatively cooled graphite skirt. The choice of graphite 
for this skirt, extending from an area ratio of 24 to 1 to the 
exit plane (area ratio 100 to 1) , rather than the more commonly 
used refractory metals, was based upon the superior nucleonic 
properties of low atomic weight graphite as opposed to those 
of high atomic weight metals. While the thermo-mechanical 
properties of graphite are quite adequate , fc or its application 
as a skirt material, there exists however a question of its 
long term compatibility with the high temperature hydrogen 
propellant. 

It is well established (Ref .3-2, -3, -4) that at high temperature 
(1000 °K) and particularly under conditions of high temperatures 
and pressures (0 (1) atm) hot graphite will react with hydrogen 
to form simple hydrocarbons. Such a reaction, if sufficiently 
vigorous, would have a pronounced effect upon the utility of 
the nozzle skirt for multiple missions requiring cumulative 
burn times on the order of 10 hours or sc. Consequently, it 
was deemed of importance to have definitive information avail- 
able concerning the details of hydrogen-carbon surface reaction 
under conditions anticipated to be present at the surface of 
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the skirt during all phases of NERVA engine operation. 

The mechanisms by which graphite surfaces react with heated 

hydrogen have been studied basically under equilibrium conditions. 

Clarke and Fox (Ref. 3- 2) employed a filament to investigate the reactio. 

of hydrogen and graphite over a temperature range 2000° K to 
o -2 

3400 K and pressures from 10 to 1 atms. These experiments, 

conducted at low total hydrogen flow rates, permitted them to 

conclude that at temperatures below 3000®K and pressures be- 
-2 

tween 10 and 1 atms, the rate of reaction at the surface of 
the graphite is directly proportional to the hydrogen pressure 
and to the square root of the dissociation constant for this 
gas. At higher temperatures the sublimation rate of graphite 
is controlling, and C 2 H 2 and C 3 H^are predominant as reaction 
products . 

Sanders (Ref. 3-2) of Lewis Research Center measured weight loss 
from graphite test samples heated in the range 1820°K to 2430°K 
at atmospheric pressure. The flow rate of hydrogen was low 

3 

(47 cm / sec) and equilibrium was maintained throughout the 

test series. The data collected from these tests at temperatures 

o 

xn excess of 2070 K, correlate quite well with thermodynamic 

calculations of the equilibrium reaction of carbon and hydrogen. 
Divergence between the data and computed reaction rates in the 
temperature range between 1820*1? and ?070°K was attributed to 
a lack of thermodynamic equilibrium induced by the low, but 
significant, hydrogen flow rate employed during the tests. 

In addition to predicting the reaction rate, the calculations 


indicated that acetylene (C 2 H 2 ) would be the primary reaction 
product. Earlier, Heddon (Ref. 3-4) determined the rate of methane 
production for graphite reacting with hydrogen at pressures 
from 10 to 10 atmospheres and temperatures from 1270°K to 
1520 P K. It was concluded that the rate determining step in 
this pressure and temperature regime consisted of hydro- 
genation of CH 2 -CI} 2 -groups at the graphite lattice coupled 
with a splitting of the -CH 2 -CH 2 - bond. It was further con- 
cluded that the reaction rate was negligible at temperatures 
below 1270®K. 

These studies, all conducted under equilibrium conditions, 
have but limited applicability to the problem of establishing 
the erosion rate of a contoured r.ozzle. It is necessary, when 
trying to establish the chemical erosion rate of the graphite 
nozzle skirt in question, to not only have knowledge of the 
graphite-hydrogen reaction mechanism, but also to have a fully 
developed model for the flow field and its interaction with the 
boundary. 

The nozzle skirt under continuous operating conditions will 
attain an equilibrium temperature which varies as a function of 
longitudinal position, assuming axial symmetry, and is determined 
by a balance between the heat flux to the surface through. the 
reacting boundary layer inside and thermal radiation from the 
outside surface. Consequently, the conditions existing at the 
graphite-hydrogen interface will be very much different from 
those maintained in the experiments referenced. In particular, 
it is difficult to determine the influence surface temperature 
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has upon the overall reaction rate from, for instance, the 
effect of hydrogen temperature or velocity. Such information 
is necessary for the proper development of a coherent model 
of the reacting boundary layer in the skirt. 

A practical engineering assault on the problem of graphite 

skirt erosion (Ref. 3- 5, -6) has shown tnae,on the basis of model tests 
employing AGCarb- 101 material as a test sample, and pressures 
and temperatures commensurate with those anticipated to be 
encountered, erosion rates will be acceptable. Similar tests 
conducted with helium rather than hydrogen failed to cause 
erosion. It might be concluded that this is an unambiguous 
indication of a purely chemical, as opposed to mechanical, 
erosion mechanism. However, the question can be raised as to 
whether the surface is not acting as a catalyst for the re- 
combination of the dissociated hydrogen present in the flow 
leading to intensified reactivity in excess of that to be 
anticipated from purely thermal considerations. A more exten- 
sive series of tests employing other gases (e.g., N 2 ) would be 
enlightening. 

While the model tests seemingly appear to provide a solid 
basis for the verification of the suitability of the existing 
skirt design, i.e., indicated acceptably low erosion rates, 
these tests have constantly failed to duplicate the low pressure 
high velocity and Mach number conditions to be experienced by 
the skirt. It is not possible at this point, on the basis 
of these or the earlier tests to predict with any assurance 
what the erosion rate will be in the actual nozzle skirt flow 
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regime. This statement is particularly true with regard to 
future higher temperature versions of the NERVA engine. As 
indicated earlier, to do so would require a comprehensive 
description of the reacting boundary layer coupled with a real- 
istic description of the reaction mechanisms present when 
hydrogen impinges on a heated graphite surface. 

It was the purpose of the experimental program undertaken 
during this contract to provide information concerning the 
detailed mechanisms by which hydrogen reacts with solid 
carbon (graphite) surfaces. 
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Experimental Program 
Approach 

In order for any experiment directed toward the invest- 
igation of the interaction of hydrogen with graphite surfaces 
to provide meaningful results it is mandatory that the basic 
physical characteristics of both the hydrogen gas and the 
graphite surface be known with precision and be independently 
controlable over a substantial range of values. An experiment 
so designed would then permit the pertinent parameters associated 
with the surface reactions to be independently varied and their 
effect upon the overall reaction schematic to be unambiguously 
assayed. 

The use of a molecular beam nozzle generating a precisely 
controlled jet of hydrogen having known physical properties 
impinging on an accurately controlled graphite surface was 
chosen as the most straight forward means to satisfy the basic 
objectives of the study. 

Experimental Apparatus 

The experiment consisted then of four basic elements: 

A. Hydrogen molecular beam nozzle 

B. Graphite surface and its holder and power supply 

C. Diagnostic instrumentation 

D. Vacuum facility in which the experiment was conducted 

The following brief description of each of these components 
serve to detail the salient characteristics of the experiment. 

• A. Hydrogen Molecular Beam Nozzle 

Figure 1 is a photograph of the hydrogen molecular beam 
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nozzle. In the nozzle design gaseous hydrogen is introduced as 
shown andjby splitting the flow^ is caused to enter both ends 
of the resistively heated rhenium nozzle tube. A small hole 
{3^m diameter) drilled midway between the ends of the tube, 
along a radius, forms the nozzle proper. 

By controlling the stagnation pressure and the stagnation 

temperature (by adjustment of the current level through the 

tube) a wide range of flow conditions could be obtained. 

Stagnation pressures of up to 20 atmospheres at stagnation 

temperatures up to 2800°K were possible, at flow rates up to 

a maximum (determined by the facility pumping capacity) of 
20 

10 part ic le s/sec . 

The Mach numbers of the beam generated by such a 
nozzle were sufficiently high ( (0)M=10) so that for the purposes 
of the experiment uhe energy distribution could within narrow 
limits (ia-10%) be considered mono-enegetic . The use of an 
almost mono-energetic beam allows the influence of particle 
velocity on the overall reaction kinetic rate to be precisely 
determined. 

• B. Grajphite Sample and Holder 

Candidate graphite samples were mounted in a sample holder 
which would permit the surface to be oriented over 2 TT steradians 
with respect to the center line of the beam. By varying the 
beam impingment angle the influence of surface anisotropic it ies 
upon the surface kinetics could be weighed. In addition 
provision was made for the samples to be resistively heated to 
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predetermined, controlled temperature levels to allow the 

influence of surface temperature to be compared with the effects 
of hydrogen particle energy. 

• Diagnostic Instrumentation 

A series of increasingly sophisticated instrumentation 
was available for this study. Initial studies of surface 
erosion under various beam and surface conditions were to be 
conducted by qualitative and quantitative evaluation of the 
amount of material removed by the reaction process. By using 
these data as a guide a detailed investigation using a mass 
spectrometer to assess the reaction products would then permit 
a detailed evaluation of the specific surface reactions 
contributing to surface erosion. 

• D. Vacuum Facility 

The use of a high vacuum environment was essential to 

this program. By maintaining the background pressure in the 
-4 

rangeflO torr two important advantages were obtained; the 
absence of significant extraneous material that could mask 
the effects of the reaction was assured and the reaction could be 
observed free from the influence of additional interactions 
due to multiple collisions in the vicinity of the surface. 

The pumping capacity associated with the vacuum facility limits 
the maxium nozzle flow rate that can be used. 

All experiments conducted in this program were run in 

the Guggenheim Physical Laboratory facility which permitted pressure 

-4 

°f "10 torr to be maintained at particle flow rates of up to 

2 o 

10 sec . Pressure levels below 10 ’ torr correspond to mean 
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free paths in excess of one meter insuring favorable conditions 
for the observation of the reaction at the surface. 

It should be noted that, in*all but one respect, the 
physical characteristics of the hydrogen beam are quite similar 
to those anticipated to be present in the nozzle skirt of the NERVA 
flight engine. The one exception to this correspondence is the 
particle density. The hydrogen density at the surface of the graphite 

test sample, when being impinged upon by the beam, was on the order 

14 3 16 17 3 

of 10 1 to 10 1 particles/cm as compared to 10 to 10 particles/cm 

in the actual nozzle. The simulation of nozzle skirt surface operat- 
ing conditions was not the goal of this work. It is purely coinci- 
dental that the experimental parameters approximate those anticipated 
to be present. 

Overall Progress 

Considerable effort was devoted, during the early phases 

of the program, toward the development of a nozzle capable of 

providing the requisite hydrogen flows at the pressures and 

temperatures required for the test program. Initial attempts 

to produce a nozzle beam using tantalum nozzle tubes revealed 

e 

that operation at modest temperatures (ie, (0) 1000 K) soon led to 
catastrophic failure because of embrittlement. As a 
consequence the use of rhenium nozzle tubes was pioneered. 

A supply of vapor deposited rhenium, and tungsten- 25%, 
rhenium nozzle tubes were secured and tested for long term 
compatibility with hydrogen under conditions of interest for 
this experiment series. The rhenium tubes supplied by The 
Advanced Rocket Technology Corporation of Irvine, California 
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i 

i 

i 

i 

! 

proved to be more than adequate to the task and superior to the 

j 

more brittle tungsten/ rhenium alloy tubes tested. j 

The use of rhenium as a nozzle material necessitated the 
development of special drilling techniques to produce the 
microscopic holes required for the generation of high Mach 
number beams. By using specially ground tungsten carbide drill 
bits and a microdrilling apparatus purchased from National 
Jet Company of Cumberland, Maryland a number of nozzles were 
fabricated for use in the program. 

An extensive series of tests were conducted to accurately 
determine the velocity distribution of the hydrogen beam 
produced by these nozzles when operated over a wide range of 
stagnation pressures and temperatures. These tests, using 
existing time of flight velocity distribution apparatus, revealed 
the Mach number to be in the range of 10 to 20 for all operating 
conditions and therefore adequate for the initial series of 
erosion tests. 

The first series of tests had as its goal the qualitiative j 

evaluation of the various molded graphite skirt materials, 
listed in Table 1. For these tests the unheated 2 cm x 4 cm 
polished samples were mounted normal to the beam center line 
and 5 cm from the nozzle exit plane. Meticulous attention 
was devoted to assuring that the test conditions were maintained 
1 the same for each sample. By operating at a stagnation temp- 

erature of 2000 °K and a pressure of 3 atmospheres hydrogen 
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particle mean energy of 0.6 ev was obtained. Each sample was 
irradiated for a period of three (3) hours, during which the flow 
particle energy was maintained constant within 2 %, 

A thin tantalum shadow strip was afixed across the center 
of the sample piece to facilitate comparison of irradiated and 
unexposed surfaces. Figure 2 is a photograph taken in visible 
light of an irradiated sample, in this photograph the light 
horizontal band is the un irradiated surface, the darker areas 
having been exposed to the hydrogen beam. The small light 
rectangular areas at the corners were shielded from the 

impinging beam by clips used to restrain the sample in the 
holder. 

The sample shown in this photograph (Sample 172) was 
visually quite similar to the other seven tested. In all 
cases visual inspection failed to reveal any marked differences 
amongst the samples studied. The demarcation between irradiated and 
unirradiated areas of the samples studied was sharp and quite 
distinct even for the areas shielded by the clips whose edges 
were raised somewhat above the sample surface, (the tantalum 
shield was in direct contact with the sample surface) . 

Because of this unambiguous visual evidence of beam 
impingement induced surface alteration it was anticipated that 
photomicrographs taken with a scanning electron microscope 
(SEM) would permit a quantitative measurement of surface 
erosion to be made. The strategm employed was to orient the 
sample such that the step in the surface (viewing from the eroded, 
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irradiated surface toward the shielded surface) could be 
clearly imaged and its elevation measured. 

Fig ire 3 is a photomicrograph of the boundary region 
between the irradiated area (right and below) and the shielded 
area (left and above) of sample 172. This material displayed 
the highest degree of surface alteration of all the samples 
studied with the SEM. In fact it was the only sample of the 
eight (8) tested in which it was possible to detect the boundary 
region at any level of magnification. 

This result was most unexpected because of the except- 
ionally sharp, distinct visual evidence all of the samples 
displayed for beam induced alteration. Inspection of figure 3 
and a similar photograph of this sample surface obtained at 
a manification of 4000X (displayed as figure 4) showed that the 
radiated surface has undergone what can best be described as 
mild excoriation. Surface irregularities in the irradiated 
area are observed to be softly rounded and generally subdued 
as compared to the shielded zone. 

Close scrutiny of the photomicrographs taken of the eight (8) 
test samples at various view angles and at magnifications up to 
4000X reveal, at most, marginal changes in the surface textures 
of all except #172. On the basis of this evidence the sample 
#172 material can be considered to be more susceptible to attack 
by energetic hydrogen atoms than the others tested. 

This conclusion must be mollified by the realization that 
the surface changes apparent could very well have resulted from 
a process that was initially quite pronounced during the three (3) 
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hour test period and that no or very little change occured 
thereafter. The resolution of this question is a moot point 
pending further tests involving both shorter and longer 
irradiation times with higher surface temperature, higher 
beam particle energies and varied sample orientations. 

Although planned these tests were not implemented during 
this contract because the overall emphasis of the program was 
shifted soon after the data were collected and in the process 
of being analyzed. At the conclusion of the reported sample 
test series major effort was devoted toward the refinement of 
the rhenium, high pressure, high temperature nozzle design. As 
a result of these refinements to the nozzle design stemming from 
an extended series of tests at temperatures in the 2600 to 2700 C K 
range at pressures in excess of one atmosphere total nozzle 

C 

longenity has been extended to 30 hours at 2700 K. 

On the basis of the known crystallographic and sublimation 
properties of rhenium a 30 hour useful life at 2700 K corresponds 
to virtually indefinite ( 1000 hours) service life at operating 
temperatures below 2400 ®K. As a consequence a high temperature, 
high pressure molecular beam nozzle is now available for the 
acceleration of all non-oxidizing gases and vapors, a signal 
advance in the field of molecular beam research. 


TABLE I 


* 


Manufacturing Historie s of Molded Gr aph ite ( 1 ) 
TEST SAMPLES 


Sample 

Number 

Filler 

Binder 

Heat 

Treatment 

£l_ 

143 (2 > 

POCO PCD (-OQ) 

cai< 4 > 

2800 

152 <2> : 

CAI (1000 C) 

ITX* 4 * 

2800 

162 * 1 2 ' 3 4 * 

ITX (1000) 

ITX 

2800 

169 (2) 

POCO PCD (-9) Q 

ITX 

2800 

171 (2) 

POCO PCD (-0) Q 

ITX 

2800 

172 

POCO PCD (-9) Q 

Ashland 
240 Pitch 

2800 

175 

ITX (700*0 

ITX 

2800 

POCOAXF5Q 

Commercial - No History 

2500 


(1) Filler and binder molded at 1600 psi and 1400 °C, except 
for POCO AXF5Q. Heated without pressure to 2800 °C. 

(2) Impregnated with CAI Polymer- 

(3) Same graphite as 175 except 162 was impregnated. 

(4) CAI - Cinnamylideneindene. 

NOTES : 

1. These materials were made for NERVA Program# SNC-95. 

2. Supplied by Mr. J. M. Napier# Nuclear Division# Union 
Carbide Corporation P.0. Box Y, Oak Ridge, Tennessee 37330. 


ETHANOL COOLING 
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ORIGINAL PAGi: k- 
OF POOR QUALITY 




T 
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Figure 3. Scanning Electron Microscope Photograph 
between irradiated surface (on right) and shielded 
(magnification X2000) ), Sample #172 


of Boundary 
surface 



t 


Figure 4. Same as figure 3 but magnification X4000 


53 


REFERENCES 


1-1 Sharma, 0. P., "Status Report on High Speed Hydrogen/Graphite 
Interaction," NGR 31-001-230, January 23, 1973, AMS Dept. 
Princeton University. 

1-2 Herring, H. J. and Mellor, G. L. , "A Computer Program to Cal- 
culate Laminar and Turbulent Boundary Layer Development in 
Compressible Flow With Heat Transfer," Princeton University, 
Dept, of Aerospace and i'echanical Sciences Report, June 1972. 


2-1 Frey, H, M. , Nickerson, G, R., and Tyson, T. J. , "ODK, One- 
Dimensional Kinetic Reference Program," Dynamic Science, 

Irvine, September 1970. 

2-2 Nickerson, G. R., "Striated Transonic Flow in a Converging- 
Diverging Nozzle," Dynamic Science, Irvine, February 1970. 

2-3 Levine, J. N., e.a., "Transonic Flow in a Converging-Diverging 
Nozzle ' Dynamic Science, Irvine, September 1970. 

2-4 Der, J. J,, "Theoretical Studies of Supersonic Two-Dimensional 
Axisymmetric Non-Equilibrium Flow, Including Calculations of 
Flow Through a Nozzle," NASA TR R-164, Ames Research Center, 
Moffett Field, December 1963. 

2-5 Frey, H, M. and Nickerson, G. R., "TDK, Two-Dimensional Kin- 
etic Reference Program," Dynamic Science, Irvine, December 1970. 

2-6 Dergazarian, T, E., e.a., "JANNAF Thermochemical Data," Thermal 
Laboratory, Dow Chemical Company, Midland, December 1960. 

2-7 Vincenti, W. G, and Krueger, C. H,, "Introduction to Physical 

Gas Dynamics," John Wiley and Sons, Inc., New York, August 1970. 

2-8 Moretti, G. and Abbett, M. , "Non-Equilibrium Effects for Quasi 
Equilibrium Flows in Nozzles," AIAA Paper No. 66-626, AIAA 
Second Propulsion Joint Specialists Conference, Colorado Springs, 
June 1966. 

2-9 Kliegel, J. R. and Tyson, T. J., "Discussion of Chemical Kinetic 
Integration Techniques," WSCI Paper No. 68-45, Combustion 
Institute and University of California, Session on Chemical 
Kinetics Calculations, Berkeley, October 1968. 

2-10 Eschenroeder , A. Q. and Lordi, J. A., "Catalysis of Recombina- 
tion in Non-Equilibrium Nozzle Flow," 9th Symposium of Combus- 
tion 

2-11 2eleznik, F. J, and Gordon, S. , "A General IBM 704 or 7090 

Computer Program for Computation of Chemical Equilibrium Compo- 
sitions, Rocket Performance and Chapman-Jouquet Detonations," 

NASA TN D-14S4; Lewis Research Center; Cleveland; October 1962. 


L 


- 54 - 

REFERENCES CON'T. 


3-1 Altseimer, J. H., Mader, G. F., and Stewart, J. J., "Operating 
Characteristics and Requirements for the NERVA Flight Engine," 
A1AA Paper No. 70-676, A1AA 6th Propulsion Joint Specialist 
Conference, San Diego, California, ,June 15-19, 1970. 

3-2 Clarke, J. T. and Fox, B. R., "Reaction of Graphite Filaments 
with Hydrogen Above 2000 K," J. of Chem. Phys., 46, 827, (1967) 

3-3 Sanders, W. A., "Rates of Hydrogen-Graphite Reaction Between 
1550 and 2260 C," Lewis Research Center, NASA TN D-2738, 

April 1965. 

3-4 Heddon, K. , "The Formation of Methane from Hydrogen and Carbon 
at High Temperatures and Pressures," Proc. of the Conf. on 
Carbon, p. 125, 1961. 

3-5 Brown, R. L. , "HES AGCarb-101 Erosion Test Heat Transfer Analy- 
sis," Aerojet Nuclear Systems Co., Sacramento, Internal Memo, 
to: Cloosen, L. B., 8 June 1970. 

3-6 Williams, J. J., "AGCarb Erosion/Corrosion Sealing," (WO 1188- 
26-524) Aerojet Memorandum, to: Shurly, L. A., 8 June 1970. 


1 


- 55 - 


APPENDIX A 


THEORETICAL DERIVATIONS 


Considering a mixture of hydrogen atoms and molecules, and 
defining as the ratio of the mass of hydrogen atoms oresent 
in the mixture to the total mass of the mixture, the number 
densities of'H and H 2 can be determined as 


• n H 


2 ot 


m. 


2m, 


and 


1- «C 1- oc 


m 


H, 


2m 


'H 


where mg is the mass of a hydrogen atom. 

Using these exoressions the equation of state can be written 

where k is the Boltzmann constant. 

Ftora equation (Al) can be deduced 


(A1) 


5- aH u “ )T - 


m.w. 


mixture 


0r 9 as 


2w 


m.w. 


H 


mixture 1 +td 


where Wg is the« atomic weight of hydrogen and 3S is the uni- 
versal gas oonstant, it follows 


* 

w. 


k_ 

m. 


(A2) 


H H 
The enthalpy is defined 
h ■ + e 

where e is the soecifio internal energy of the mixture, with 
•ero energy taken to be that of the atomic species at absolute 
sero temperature. 

As only vibrational, rotational, translational and dissociation 
energies of the speoies are considered (electronic excitation 
is negleoted), e can be written 
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2l n i ^dissociation^ translation + ^vibration^ 

x g \ 

rotation^' 

summed over both softoies i in the mixture. 

Translation, rotation and vibration are assumed to be in 
equilibrium, i.e. they have a Maxwell-Boltzmann distribution. 

Expressions for the energies are then 

B rot * “ 

E 1 ke Vl 

vi^" ^ 7 ? 
e 

* tP - 3/2 kT 

E dis i * ke Di 


where 


®Vi * characteristic vibrational teraoerature of 


species i 


e Di * characteristic dissociation temperature of 
soecies i. 

As the atomic soecies has only translational energy, one 
obtains for the internal energy 

• ’ jib ( * 5/2 * -Sf2S-( V '(. V % * e B ))T (A3) 


and for the enthaloy 


h • 2^ (3 “ /i! * 7/2 ♦ i f !S (e v /(e 8,r/T - 1 ) ♦ e B ))T (A4) 
or 

e - e(< ,T) 

h « h(«4,T) 

®V and 0 D are the characteristic 1 temperatures of molecular 
hydrogen. 

For the calculation of the soeoific heat ratio y and of 
some other quantities we will need the partial derivatives 
of the enthaloy. 
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where 


«*•! 


'He*''!,) 


. * fiV 
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(A*> 

(AO 


Ae the specific heat at constant volume Cy is defined 


«.„= 2 . o . 2 ) 


r ^ t - R 


where 




lUt.t 


is the specific gas constant of the mixture, and 


v: br « -i + iL 

«V Cy 


it follows 


r-- 


jc ClifO 

a*n M u.,. 


For the local frozen a peed of Bound can now he written 


ft} * \fy RT a { Ji y [ 

1vm M 1 


t A 7 .) 


(AS) 


He combi nation linte 


Considering the two third body reactions 
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where 


k f * forward reaction rate of the i th third body reaction 
backward reaction rate of the i tb third body reaction, 

the dissociation rate becomes (c.f, ref, 2-7) 




gi 1 ) 




where 


K c ■ JT” ■ the equilibrium constant based on concentrations per 
i unit volume for the reaction H,,~ 31. 


As the third body has no effect on the equilibrium concentrations, 

write 

k # ■ K „ K and k- - K k. 
f i c b l *2 c \ 

Substitution of these expressions in equation (A9) yields 


one can 


16* > i Ct'-O - ^ 1 

W H ICt. 4 


(A to) 


Ref. 2-10 gives for the values of and ^ the empirical expressions 
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Substituting in the equation f or to , one obtains 


% ^ L V„K t J 


^[l'— J*« - 22S?] 


As ref. 2-6 gives the equilibrium constant based on partial pressures, we will 
have to convert the last equation with the aid of 


lc- . -he? . J^2_heh 


3?T L-T i 


where is the equilibrium constant based on partial pressures, which results 


, £hj- 

iw H ' 


-v-') 


Eliminating^ with the equation of state (Al) it follows 

k.Tt» l* <|p J 


For equilibrium flow and chemically frozen flow the rate of dissociation 
is equal to zero, which gives us for frozen flow 

f ■ constant 


and for equilibrium flow 
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hence 


■a « ^ 


ilE 

S‘ 


CaiO 


As is a function of T only, is a function of p and T only. The 
partial derivatives with respect to p and T are 
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One-Dimensional Non-Equilibrium Flow 


For one-dimensional flow the equations of conservation of momentum and 
energy and the chemical relaxation equation can he written 

u_ £T v b* « . i «*£ „ 0 

^ = o 

do c. doc. 


cUc 

dot 


to 

w 


as shown in 2. 1.1. 2. 


A fourth differential equation for p can he derived by differentiating 
the equation of state (Al) 

•5 , r ‘ dp t dT | cU I 0 
«l% ^ P dx ’ T dx "iMtcTocl^ 


Rearranging these equations 
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Combining equations (Al) and (A7) X can be written 
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Substituting this result in the expressions for the functions f , , f . f 

1 * 2^ 1 

and f. and rearraneinc them 


\ and rearranging them 
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For the computation of the one— dimensional non-equilibrium flow we will need 
the partial derivatives of the functions f^, B and A with respect to x, w, 

*< , T and 

Defining 
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and setting 
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these derivatives can easily be shown to he 
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Ih« derivatives y T and ^ of the specific heat ratio can he shown to be 
from equation (A?) 
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